Bang-Bang Control Design for Quantum State 
Transfer based on Hyperspherical Coordinates and 
Optimal Time-energy Control 

Weiwei Zhou 1 , S. G. Schirmer 2 , Ming Zhang 1 , Hong-Yi 
Dai 3 

department of Automatic Control, College of Mechatronics and Automation, 
National University of Defense Technology, Changsha 410073, People's 
Republicof China 

2 Dcpartment of Applied Mathematics and Theoretical Physics, University of 
Cambridge, Cambridge, CB3 OWA, UK 

3 Department of Physics, College of Science, National University of Defense 
Technology, Changsha 410073, People's Republic of China 

E-mail: zwwannSnudt.edu.cn, sgs290cam.ac.uk, zhangming@nudt.edu.cn, 
hydai@nudt . edu . cn 

Abstract. We present a constructive control scheme for solving quantum state 
engineering problems based on a parametrization of the state vector in terms 
of complex hyperspherical coordinates. Unlike many control schemes based on 
factorization of unitary operators, the scheme gives explicit expressions for all 
generalized Eulcr angles in terms of the hyperspherical coordinates of the initial 
and final states. The factorization, when applicable, has added benefits that phase 
rotations can be combined and performed concurrently. The control procedure 
can be realized using simple bang-bang or square- wave-function controls. Optimal 
time-energy control is considered to find the optimal control amplitude. The 
extension of the scheme to implement arbitrary unitary operators is also discussed. 
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1. Introduction 

Control of phenomena governed by the laws of quantum mechanics is increasingly 
recognized as a crucial task and prerequiste to realizing promising new technologies 
based on quantum effects from the use of photonic reagents in chemistry pQ to quantum 
metrology and quantum information processing [2] to mention only a few examples. 
From early beginnings in the 1980s (see e.g., [3J 0] |SJ HI [7] ) , there has been considerable 
recent progress in both theory and experiment of quantum control [3] ■ Among the 
core tasks for quantum control are quantum state and operator engineering. In the 
former case the main objective is to prepare the system in a desired state, which is 
usually a pure state \ijjf) represented by a unit vector in a Hilbert space H associated 
with the system. The task can take various forms, from state transfer, i.e., steering 
the system from a known initial state |^o) to the target state [TUl EJ H21 E3 HI], 
to purification or state reduction, i.e., preparation of a desired pure state starting 
with a mixed or unknown initial state, usually involving some form of feedback from 
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measurements of an observable p~5j [TH [T7], to protecting or stabilizing a desired 
state in the presence of environment noise or disturbance [T5J [TU [20]. Often the 
goal is the preparation of a non-classical state such as a Greenberger-Horne-Zeilinger 
(GHZ) state [5T], or a maximally entangled Bell state [55]. Operator engineering 
usually involves engineering the dynamical evolution to realize a particular unitary 
operator [23, 24, 25, 26, 27 , and plays a crucial role in the implementation of quantum 
gates in the context for quantum information processing. 

Some of the tasks mentioned above such as purification or stabilization generally 
require measurements and feedback [TS] HE] [T7] , or possibly coherent feedback [3TJ [35] , 
and some proposed control strategies actively take advantage of environmental 
effects [28] [29] [30] or even backaction effects of measurements and feedback [18] [19] . 
Perhaps the majority of control strategies for quantum dynamics, at least to date, 
however, rely on coherent open-loop control, i.e., manipulation of the dynamics via 
coherent interaction of the system with external fields or potentials, the type of 
control that we will be considered in this article. The main reason for foregoing 
measurements and feedback is to avoid the disturbance of the system that results 
from the backaction of measurement and feedback on the system, leading to complex 
non- unitary dynamics and decoherence. The challenge of open- loop control is to design 
external fields or potentials acting as controls off-line based on a model of the system. 
The main strategies for control design in this context are based either on geometric 
ideas, or more formally Lie group decompositions, as in [5J2 [511 [55] [5B] [57] , so-called 
model-based feedback design [331 EH 023 [36l [37j [38] , or formulating the problem as 
an optimal control problem and using optimization techniques (see [39] and references 
therein, [HI 05]). The latter approach has been used successfully to find solutions 
for many different types of control problems — from control of vibrational modes via 
ultrashort laser pulses 40] , to control of nuclear spin systems 0TJ 05J , to control of 
spatially distributed systems in 43 , to implemention unitary operators [44] [45] and 
encoded logic gates [46], to optimizing state transfer in spin networks [47], to the 
creation of various types of entanglement (5TJ [55] . Optimal control is important and 
holds considerable promise of enabling robust control of complex, imperfect systems 
with limited control 03] 051 05] . Often it requires control with complex temporal and 
spectral profiles, however, which may be difficult to implement for certain systems, e.g., 
in solid-state quantum dot systems controlled by voltages applied to gate electrodes, 
where it may be difficult to implement complicated time- varying voltage profiles. 

For these reasons constructive control schemes that require only simple pulses such 
as approximately piecewise-constant functions (bang-bang controls) remain useful as 
an alternative, which can often be optimized to mitigate limitations of the system or 
the control to some extent, as in the case of optimized Euler angles to compensate 
for non-orthogonal rotation axes [481 E2]i f° r example. This is the type of control 
considered here. Specifically, we consider state-transfer tasks in which initial and 
finial states are given. We show that parametrization of the initial and target states 
in hyperspherical coordinates |50j yields a simple constructive control scheme for 
state-transfer tasks that requires no complex calculations of the control parameters, 
i.e., all control parameters are given in terms of simple functions of the initial and 
final state coordinates. The scheme has some additional advantages over alternative 
geometric schemes, e.g., based on decomposition into Givens rotations [53] in that 
many operations can be performed either sequentially or in parallel, reducing the 
time required to implement the control schemes. We introduce a parameter A which 
represents the ratio of costs of time and energy, and further explore the trade- 
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off between time and energy optimal control using time-energy performance index 
J = f Q f [A + E(t)]dt, where E(t) is energy cost of Bang-Bang control at t, tf is 
terminal time. It is shown that the product of the terminal time t*j and the energy 
cost E* for optimal bounded or unbounded piecewise constant controls only depends 
on the geometric parameters of the initial and target states and is independent of A 
but A determines the optimal field strength of the controls, = V\- The scheme 
can be generalized to implement arbitrary unitary operators, and we again find that 
the resulting decomposition has some advantages in that many operations commute 
and can be performed in parallel. 

2. Pure-state Transfer by Bang-Bang Control 

Pure-states \ip) of a quantum system defined on a complex Hilbert space H with 
dimH. = N < oo can be represented by complex vectors c e C N by choosing a 
suitable basis {\n}}^ =1 for H, 

N 

|V> = ^>». (1) 
n=l 

The modulus squared |c„| 2 of the coordinates can be interpreted in terms of 
probabilities provided c is a unit vector. For most applications the global phase of the 
state is irrelevant, i.e. we can further identify \i[>) ~ e^lip). Given these considerations, 
physically distinguishable pure states can be uniquely identified with elements in the 
complex projective space CP^ -1 = S 2N ~ 1 /Si, and we can uniquely represent pure 
states by unit vectors in C N if we fix the complex phase of one coordinate. 

Pure-state transfer is the task of transforming a given pure quantum state \tp^) 
to a desired pure quantum state \tp^) and is one of the most fundamental tasks 
in control of quantum systems. Many of the control strategies mentioned in the 
introduction have been applied to this problem, including constructive control schemes 
based on the Lie group decomposition. Indeed, it is quite straightforward to see how 
to solve the state transfer problem for an TV-level system in principle, if we are able to 
implement unitary gates on a sequence of connected two- level subspaces. Assume, e.g., 
that SU(2) operations can be implemented on the subspaces spanned by {|1), |2)}, 
{|2), |3)}, . . . , {\N — 1), \N)}. We can decompose any unitary operator in SU(JV) into 
a sequence of SU(2) rotations on these two-dimensional (2D) subspaces. Each of these 
can be further decomposed into a sequence of three rotations about two orthogonal 
axes using the Euler decomposition. It therefore suffices if we can implement rotations 
about two fixed orthogonal axes on each of the 2D subspaces. Applied to the problem 
of quantum state transfer, it is not difficult to see that we can transform any complex 
unit vector c^°) into any other complex unit vector ^ by a sequence of N— 1 rotations 
on the 2D subspaces defined above 

C< s ) = [/( JV_1 ' jV ) . . . [7(2,3)^(1,2)^0) (2) 

where ;/( n >™+ 1 ) indicates a complex rotation on the subspace spanned by {|n), |n + l)}. 
Decomposing each U^ n ' n+1 ^ further into three rotations about two fixed orthogonal 
axes, u[ n ' n+1 \a) and U^^^ift), by suitable angles jk, 



jjKn+i) = ui n ' n+1) (y 3 )U^ n+1) (^)U^' n+1) ('n), 



(3) 
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shows that in general 3(N — 1) such rotations are required to transform a given initial 
state to a target state using a sequence of elementary unitary transformations, 

XU^MU^\ 72 )U^M^. (4) 

It is easy to see how to transform pure states in principle, but it is not obvious how 
to derive the correct rotation angles "fk in the sequence, which is what matters in 
practice. Although it is possible to constructively compute the jk, the dependence of 
7/c on the state vectors c^°) and c^^ is complicated. 



3. Bang-Bang Control Scheme based on Hyperspherical Parametrization 

In this section we discuss how to obtain explicit expressions for the rotation angles jk 
and show that it can be easily solved by parameterizing the initial and target states 
in terms of complex hyperspherical coordinates. 



3.1. Complex hyperspherical coordinates 

Any complex unit vector c can be parametrized in terms of complex hyperspherical 
coordinates (9, 6), 



( C1 ^ 




/ COS#i 


C2 




sin6»i cos 9 2 




= e^ 




CN-1 




e i0N-2 s - m Q 1 _ . . s in6» A r_ 2 C0s6' A r 


\CN ) 




\ e^™- 1 sin 6>i . . . sin6>jv-i 



(5) 



where 9 and 4> are vectors in with < 9 n < f and — ir < (j) n < ir, and e 1 ^ 

is a global phase factor, which is usually negligible. Thus, assuming normalization 
and neglecting global phases, any pure state is uniquely determined by its complex 
hyperspherical coordinates (9, $) which can be calculated easily by Algorithm [lj 

Although there are many equivalent parameterizations of pure state vectors, the 
beauty of complex hyperspherical coordinates is that we can easily give an explicit 
constructive bang-bang control scheme for state transfer l^ )) i— > \ip^) such that all 
control pulses are determined directly by the coordinates of the initial and final states 



3.2. Control Assumptions 

The following scheme is based on the assumptions that (a) we can neglect free evolution 
Hq = 0; (b) we have local phase control, i.e., we can implement control operators that 
introduce a local phase shift, 

Z n = U n , n = 2,...,N (6) 

where 1^ is the identity on % and II„ is the projector onto the subspace of % spanned 
by the basis state |n); and (c) we can individually control transitions between adjacent 
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(6, (f>) <- HyperCoord (c) 
Compute complex hyperspherical coordinates 
In: c complex vector/pure state 
Out: 6, cf> hyper-spherical coordinates 



1 


N <- length(c) 


2 


c <- c/norm(c) 


3 


c «- exp(— i * angle(ci)) 


4 


angle(c 2: at) 


5 


a <- ABS(c) 


6 


#i <— arccos(ai) 


7 


si <— sin(0i) 


8 


for n <- 2,...,N- 1 


9 


0„ «- arccos(a„/s„_i) 


10 


Sn <- s n _i sin(6>„) 



Algorithm 1: Computation of complex hyperspherical coordinates 



energy levels, i.e. that we can realize control Hamiltonians of the form X n or Y n , 

X n = {\n + l)(n\ + \n)(n + l\), n = 1, . . . , N - 1. (7a) 
Y n = i(\n+l){n\-\n)(n+l\), n = 1, . . . ,N - 1. (7b) 

The evolution of the system under any Hamiltonian H is governed by the Schrodinger 
equation 

ihU{t) = HU(t), U(0) = I N , (8) 

and we choose units such that the Planck constant K = 1. This shows that the 
evolution under the control Hamiltonian H 6 {LX n , LY n , LZ n } is given by the one- 
parameter groups exp(— iLtX n ), exp{—iLtY n ) and exp(— iLtZ n ), respectively. The 
evolution is unitary as the operators X n , Z n and Y n are Hermitian. In particular, this 
means that we can implement the complex rotations 

u n( a ) = exp(-iaX n ), JjJ L (a) = exp(-iaY n ), C/jf (a) = exp(-iaZ n ), (9) 

by applying the control Hamiltonians LX n , LY n and LZ n respectively for some time 
t — a/L. In the following only two types of control operations {X n , Z n } or {Y n , Z n } 
are required. 

The assumptions on the control Hamiltonian are somewhat demanding, although 
no more so than the control requirements for the standard geometric decomposition 
Eq.Q. While these requirements cannot always be satisfied, there are systems for 
which these control operations are quite natural such as a charged particle trapped in 
a multi-well potential created and controlled by surface control electrodes as shown in 
Fig. |Tj A physical realization of such a system could be a multi-well potential created in 
a 2D electron gas in a semiconductor material by surface control electrodes. Changing 
the voltages applied to different control electrode enables us to vary the depth of 
individual wells as well as the height of the potential barrier between adjacent wells 
and thus the tunnelling rate, giving raise Z n and Y n rotations, respectively. 



Bang-Bang Control Design for Quantum State Transfer 6 




Potential well Tunnel barrier 

Figure 1. Charged particle trapped in a multi-well potential created by control 
electrodes. Red electrodes allow control of potential barriers and thus tunnelling 
rates, while blue electrodes control depths of the wells 1 to 3 and thus their energy 
levels. We can choose default voltage settings such that all wells have the same 
depth and there is no tunnelling, so that we effectively have Hq = 0. Then by 
raising or lowering the voltage of electrode 1 we can introduce a relative phase 
shift between the ground state |1) of the first well and the other two ground states, 
and by changing the voltage applied to the red electrode between wells 1 and 2 
we can induce tunnelling between the first two wells, etc. 



3.3. Explicit Control Sequence 

To illustrate the constructive procedure, let us consider the case N = 3 with Y, Z 
controls. In this case the control operators take the explicit form 



Z 2 = 




Z, 




and the corresponding evolution operators are 



Uf(a) = ( 
,0 



Ui{a) = 




Y 9 




' cos a — sin a > 
Ui (a) = sin a cos a 
1, 



U?(a) 




Given these control operators and the hyperspherical coordinate representation of the 
initial and target states, it is now very easy to see how to steer an arbitrary initial 
state to an arbitrary target state in the following seven steps: 



Step 1. 



? (0) 



?(0). 
7 2 



1 










1 










e -i4 



Step 2. 



?i 0) , 



4 0) ;. 




, 0): Apply phase rotation U% 



cos 



e^i sin (9 UJ cos 9. 



? (0) 
(0), 



3(0) 



V 



9f 




0,0): Apply phase rotation U 2 
( 



COSC?^ 

M 0) sin^ 0) cosflf 



V 



cos^ 0) 
sin 6^ cos 9^ 
sinff) sinflf 



4 0) ) 
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Step 3. 




0;0,0): Apply population rotation U% (— 2 



( 




Step 4. (6^ 0) , 0; 0, 0) ->■ (9[ s \ , 0; 0, 0): Apply population rotation Uf(9[ s) - 



t 

V 

Step 5. 



cos(6>i 



sin 



?i s) ,0;0,0) 




cos 6* 



sm 




cos 6* 



sm i 



l 



?^ s \ ^2^; 0) 0) : Apply population rotation U 2 Y { 



(0)j 





sin 



cos 



0« 

cos 



}(«) 
7 2 

(a) 



Step 6. (e{ s) ,e ( 2 s) ;o,o) 



3 (»). 

7 2 i 



44 



0): Apply phase rotation U 2 {—(pi ) 



Step 7. 









V 



e**i ' sin 
sin 



cos e[ s) 

^ CO8 ( 2 S) 

[ s) sin0 2 s) 



( 



\ 



cos 9{ 
iW 1 



(<0 



sm 



cos 



sin 0, 



cos 



7 2 



Apply phase rotation U$ (—<j)^) 



The generalization to N > 3 is straightforward, as shown in Algorithm [2] Given 
a Hamiltonian of the form 

2AT-1 



H 



E 



(10) 



where H 2n -i = Z n+ i, H 2n — Y n and u m (t) are controls (e.g. voltages), the bang- 
bang control sequence given by Algorithm [2] can be implemented by applying 47V — 5 
control pulses. At the fcth step we apply a constant control field Ustk) = Lk for time 
tk = Jk/Lk, while all other controls are set to (or the voltages are set to their default 
values). Notice that in practice we cannot apply fields for negative times, thus the 
sign of -Lfc must match that of 7^.. However, if 7^ is negative and > 0, we can also 
apply a field fs(k) = ^fc f° r time tk = (jk + ^)/Lk as jk + 27r > and effects the 
same rotation. 

If X n control Hamiltonians are used instead of Y n control Hamiltonians, error 
i n_1 is created in the nth coordinate by the population rotations. So the algorithm 
need be slightly modified to correct phase factors of \n)(n = 2, • • • , N). We can achieve 



this by adding Tj(n mod 4) to the phase angles 
and the phase factor of the nth coordinate is e 1 ^" 



noting that e 



Z7r/2(n mod 4) 
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Besides giving explicit expressions for the rotation angles in the decomposition, 
the scheme has an additional advantage compared the the standard decomposition 
Eq.Q considered earlier: While the rotations in the standard factorization do not 
commute, the first N — 1 and final N — 1 phase rotations in the decomposition based 
on complex hyperspherical coordinates are represented by diagonal matrices which 
commute. This means that these operations can be applied concurrently rather than 
sequentially, leading to a potentially considerable reduction in the total length of the 
control sequence. 

(5,7) <- StateTransfer (c(°',cW) 
Compute sequence of rotations required for state transfer 
In: c(°),cj( s ) initial and target state vectors 
Out: S, 7 Bang-bang control sequence 

1: (6>(°),^°)) <r- HyperCoord(c^) 
2: (6K S ),0( S )) <r- HyperCoord(c( s )) 
3: for n <- N- 1, 

4: Append S by In - 1, 7 by ffi // Apply Phase Rotation U^ +l {<f>^) 
5: for n <- N- 1, ...,2 

6: Append S by 2n, 7 by -6»1 0) // Apply Population Rotation Uj t (-6^ 0) ) 

7: Append 5 by 2, 7 by 6^ s) - 6^° } // Apply Population Rotation UY{0[ s) - 6^ 0) ) 
8: for n <- 2,...,N- 1 

9: Append S by 2n, 7 by 9 { n s) // Apply Population Rotation U%(e { n s) ) 
10: for n <- 1, . . . ,N — 1 

11: Append S by 2n - 1, 7 by -<^ s) // Apply Phase Rotation U^ +1 (-(f) { n s) ) 

Algorithm 2: Control Scheme to achieve state transfer c<°) i-> c< s ) in 4A - 5 steps 
using bang-bang control, based on hyperspherical coordinate parametrization. S and 
7 are vectors of length 4 A — 5, whose elements are integer labels indicating the control 
Hamiltonian (m = 1, . . . , 2 A — 2) and rotation angle 7^, respectively. 

3.4- Application: Creating a multi-partite- entangled W -state 

As a simple application of the scheme, suppose first that we have N sites and starting 
with only site 1 populated, i.e., in state |1), we would like to prepare an equal 
superposition of all N sites \n) for n = 1, . . . , N: 



All we need to do is compute the hyperspherical coordinates of \tp), e.g., for A = 10 
9 = (1.2490, 1.2310, 1.2094, 1.1832, 1.1503, 1.1071, 1.0472, 0.9553, 0.7854) 



to the initial state 1 1) . This results in the following sequence of states being created 



JV 




n=l 
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n = 1 


n = 2 


n = 3 


n = 4 


n = 5 


n = 6 


n = 7 


n = 8 


n = 9 


1.0000 


0.3162 


0.3162 


0.3162 


0.3162 


0.3162 


0.3162 


0.3162 


0.3162 


0.3162 





0.9487 


0.3162 


0.3162 


0.3162 


0.3162 


0.3162 


0.3162 


0.3162 


0.3162 








0.8944 


0.3162 


0.3162 


0.3162 


0.3162 


0.3162 


0.3162 


0.3162 











0.8367 


0.3162 


0.3162 


0.3162 


0.3162 


0.3162 


0.3162 














0.7746 


0.3162 


0.3162 


0.3162 


0.3162 


0.3162 

















0.7071 


0.3162 


0.3162 


0.3162 


0.3162 




















0.6325 


0.3162 


0.3162 


0.3162 























0.5477 


0.3162 


0.3162 


























0.4472 


0.3162 





























0.3162 



Each Y n extends the superposition by one site until we are left with the desired state 
after the final step. Creating such a superposition state may not seem very interesting 
but it has an interesting application in the area of entanglement creation, for instance. 

The Hamiltonians for many systems such as interacting quantum dots or 
coupled Josephson junctions, for example, can be described to a reasonably good 
approximation by an XXZ-spin network model 



i{°~ m V n + (T m <7 n + Kcr m cr «)) 



(11) 



where is an TV-fold tensor product for which the nth factor is A and all others are 
the identity, and X, Y and Z are the usual 2x2 Pauli matrices. For a chain with 
nearest-coupling we have 7 m „ = except when n — m + 1 and for k — we have 



the so-called XX-coupling model. The Hamiltonian ( 11 ) commutes with the total spin 
operator S = Z n and decomposes into excitation subspaces for any choice of a n 
and jmn- It is easy to see that a state such as the TV-state 



v^ = ^(im--4> + im--4> + -.. + iw--4t>) 



(12) 



belongs to the single excitation subspace, as does the state |tl ■ ■ • I)- On this subspace 
the Hamiltonian (111 can be simplified. For a chain with nearest-neighbour coupling 
and k 



we obtain (up to multiples of the identity): 

H\ = y -u n z n + j n x ni 



(13) 



i.e., the Hamiltonian is exactly of the form required for our scheme. Treating a n and 
7„ as control parameters, we can use the scheme above to create a W state starting 
from the product state |1) = |tl • • ■ I)- Since we can only implement X„-rotations, we 
will need to apply the phase corrections 



N 



71=2 



-f(n-l mod 4), 



which can be applied concurrently in the final step. 



4. Optimal piecewise-constant Control and Time-energy Performance 

The bang-bang control sequence given by Algorithm [2] leaves us considerable freedom 
of choice for the controls. Choosing large control amplitudes will result in short pulse 
durations, thus optimizing the transfer time tf. However, large control amplitudes 
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may not be feasible and have undesirable side effects in terms of transfering too much 
energy to the system. We can try to optimize the field amplitude by stipulating that 
the state transfer is to be achieved while minimizing a time-energy performance index 



J = 



2N-2 



[A+ Wm{t)\ 2 ]dt 



(14) 



m— 1 



where A is the ratio factor of the costs of time and energy and A > 0. Larger values 
of A indicate a stronger emphasis on time-cost, while smaller values of A give more 
weight to the energy cost of the controls. 

If the controls can take values f m (t) £ {0, ±L} and the pulses are applied strictly 
sequentially, then the total length tf of the control sequence is 



1 
L 



N-l 



N-l 



Ei^i + i^i + E^f 



(s) 



9 ( 0)l 



n=l 



n=2 



< i [ 2( iv-iK + 2 (iv-2) 2 + 2 ]=4 r l 



(15) 



because of < 9 n < | and < \<j) n \ < it. Noting that a 2 + b 2 > 2ab, with equality 



exactly if a = b, we have 

K 



K 



J = V (A + L\)t k < V 2\[XL k t k < 2y/\t f maxL k 



(16) 



k=l 



k=l 



with equality if and only if L k = \/A- This shows that the optimal choice of the field 
amplitudes is L k = a/A, for which we have 



t* f < ^^J^ , J, = min J = 2Xt* f < \/X(6N - 7)tt 



(17) 



and the corresponding optimal energy cost is E* = J* — Xt J < ^VX(6N — 7)ir. As 
expected, as A goes to 0, t*j becomes infinite and E* goes to 0, but their product 
remains constant 



t*-E* = 



(18) 

and depends only on the geometric parameters of the initial state and target states. 

If first and last N — 1 phase rotations are applied concurrently the transfer time 
is reduced 



i 

< - 

- L 



N-l 



max 



2ti + 2(N 



max|^)| + E(^) + ^)) + |^-e 

n=2 

(2A^ + 3)7T 



(0) 



(19) 



2) - + - 
'2 2 



2L 



Setting 4>^l x = max„ \(f>^\ and 0, ( nL = 

and tx — 0mL/i, and thus we must choose L n > (fin^ /h and L n = (j>n^/tK, 
respectively for the control amplitude of the first and last TV — 1 concurrent pulses to be 



c n |</>n^|, shows that we have t x = ^max/i 
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able to implement all N — 1 phase rotations concurrently in time t\ or tx , respectively. 
Furthermore the performance index changes 

N-l K 

J < 2t' f V\ max L k + V hi + V Ll, (20) 

~ ; N<k<K+l-N ^— ' ^ k K ' 

- ~ k=l k=K-N+2 

which suggests that we can improve the performance index and reduce the energy 
cost by choosing the amplitudes of the first and last N — 1 concurrent pulses to be 
as small as possible, i.e., L n = (fin /t\ and L n = (fin /tK, an d Lk — V\ for all other 
amplitudes. 



5. Implementation of Unitary Operators 

5.1. Complex hyper spherical representation of unitary operators 
Any N— dimensional unitary operator can be represented as follows: 

N 

cr=Xy w i«iX«ji 



(21) 



where {lit.,)} constructs an orthonormal basis set in N— dimensional Hilbert space. 
From (21 1, we can find that global phase factors of all \uj)s do not affect U, so they 
can be neglected. Assuming (l\uj) is real and positive for each j, by ([5| the complex 
hyperspherical parametrization for {\uj) : j = 1, • • • , N} can be given by 



f\ui)\ 

'u 2 ) 



IjV-2 







e« 



Ii o 

o e^- 1 ) 



e (JV) 



\\un)J 



|2) 



where 



cw = (^ fc )),(^w,...,( c irtw),(^) 



(22) 



(23) 



and (cf ( fe )) T is the transpose of the vector cf W and 



Cl 



(fc) 



3 (fc) 



e Vl sin 



' "i 



cos 



,(fc) 



(AO 



e^-nt^siii 



s(fc) 



(24) 



C2 



(fc) 



-e' 



x(*0 

Pfc - 2 COS 



(AO 



sin^ 
cos 0\ ' cos 6*2 



\ -e^fc-i C os 



^^m^cos^ 



(25) 
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and so forth until 



CO 



,(fe) 



fe-i/ 



(26) 



4 (fc) = 



(27) 



for 2 < fc < N. Note that (C^Y = I k 



5.2. Realization of unitary operators 

Suppose that Y n and Z n controls as defined in eqs. ^ and (7b) are permitted for 
1 < n < N — 1 and 1 < n < N, respectively. Then U can be realized by a sequence of 
bang-bang controls 



where 



N 



T=Y[U n , U n = l[U, 



N-n+j( 



*S B) ) 



(28) 



(29) 



The Z-phase rotations (underlined) commute and can be applied concurrently. Thus 
U n can be implemented in n steps and T in N(N+l)/2 — 1 steps and the entire process 
in N(N + 1) — 1 steps. If we took the more contentional approach of factoring an 
operator U G SU(7V) into a sequence of N(N + l)/2 rotations on two-level subspaces, 
e.g., spanned by {\n), \n + 1)}, and further decomposed each of these SU(2) rotations 
into three elementary Y n and Z n rotations using the Euler decomposition, we would 
require 3iV(7V + l)/2 steps instead, and since Z n and Y n operations do not commute, 
these could not be implemented concurrently. 

The proof of the result is constructive. 
(1) Effect of each U n . Let \e^) be an arbitrary state in the space spanned by 
{|l),|2),---,|JV-n)}and 



/le (n) A 

\ e N-n+ll \ 



'7V-n+2/ 



V I 



(\N-n + l)\ 

\N -n + 2) 



V \m J 
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where is as defined in Eq.(23). That is, 



'N—n+l/ 



■■ cos 6^ \N - n + 1) + e^i ' sin cos 0^ | JV -n + 2) 



■e i *»-'sin0 1 (n) 



■ sin 



COS 



(n) 
^N-n+2 



) = sin0< n) |JV-n+l) -e 1 ^" 



■ sin 



ft a* ft w 



3 (n) 



cos 0^,™ ) |JV-n + 2) 



e^-'oos^" 5 sin^ n) 



J4> 



1 cos 



g (n) 



Sill 



3 (n) 



■ sm 



sine^inC-M^) 



3 (n) 



COS 



Jn) 



i\N-l) 



{n) 



(n) 



-N-l/ 



3 (n) 



| e W) = e <- 2 sinC-il^ - 1) - e 40 "- 1 cos^il^)- 

J7„ leaves any state |e^) in the subspace spanned by {|1), |2), ■ • • , |JV — n)} 
invariant, i.e., U n \e^) = \e^) as U n is the identity on this subspace. Furthermore, 
Sec. 3.3 shows that applying U n to \en— n +x) maps it to the basis state \N — n + 1), 
and we can verify by direct computation 



\N-2) 
□(«) 



cos 0^ 2 cos ft \N-1)- cos ft sin ^ | JV) 



,(») 



3 (n) 



3 (n) 



UJe 



(n) 



= -\N -n + 2) 



»Piv-„+2> = ^-n+iPHOM'V - n + 1) - cos^r |JV - n + 2)) 



n-2 



U \e {n) ) 
u n\ c N-l/ 



= H ^-„ +j (^ n) )(sinft|JV - 2) - cos ft |JV - 1)) 



■|JV-1) 



Unletf) = II ^-n+i^^sm^JJV - 1) - cosftliV - 2)) 

3=1 

= -|JV) 

U n |eW) = |e< n >) |e(™)) € Span{|l), • • • , |JV - n)} 

U n \e%l n+1 ) = \N -n + l) (30) 



UJe 



(n)\ 



-li) 



N-n + 2<j<N. 



(2) Effect of T. Usin g th e previous result we now show that T\u n ) = ( — l)^ n 1 '|n) 
with |u„) as defined in (22). Let = (0, • • • , 0, (ct^) T ) be a row vector of length 



N where the coefficient vector is as in Eq. (|24| and the number of zeros is N — k. 
Eqs J22}-((29| and Q give 



/|1> 



T|ui) = Ux->>U N . 



Ux->> U N ^U N \e[ N) ) =U X --- U N -x\l) = |1) 
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T\u 2 ) = U 1 ---U N a^C^ 

( I 1 ) \ 

-|2) 

Furthermore, for 3 < n < N — lwe have 
U N \u n )=a< N -^ YL 

k=N-n+2 



(\e[ N) )\ 



= -U 1 ---U N - 1 \e% f - 1) ) = -\2) 



-N-k 



C^) 



N-2 



=s (jV-n+l) 



n 1 



fc=AT-n+2 



AT-3 



LjV-fe 



C^-i^v|«n> =a^-" +1 ) I] 

fc=N-n+2 



and continuing we obtain for 3 < n < N — 1 



/ I 1 ) \ 

-|2> 



/ 11) \ 



/ II) \ 

-|2> 

(-D 2 i e r 2) > 



V(-i) 9 |eiT a V 



f/iV-„ + 3---C/jv|M„) =a( Ar -"+ 1 ) 



/ (-1)°|1> \ 

(_l)(«-3)| n _ 2 ) 
( _ 1)( „_ 2)|e (iv r+2 ) ) 

V(-i) (n - 2) i e r ,i+2 v 



• • • U N \u n ) =(-l)("- 1 )|e( l JV -"+ 1 )) 

T|u„) = U N - n+1 ---U N \u n ) =(-l)("- 1 )|n), 



Finally, we have 



JV-2 



fe=3 



c( fc ) 



/ (-l)°ll> \ 
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and thus 



T\u N ) =U 2 (0, 



-,0,(cW) 



{-1)( N ~V\N-1) 

V (-iy N -^\N) J 



U 2 (-ir^\e%) 



Thus we finally have 



N 



N N 



TUT^ =Y,e ilfn T\u n )(u n \T^ 



^e^|n)(n| = n^(^) 



n—1 n—1 



and 1/ = rt n^=i U%{<p n )T as claimed. 



6. Discussions and Conclusion 

We have presented an explicit geometric control scheme for quantum state transfer 
problems based on a paramctrization of the pure state vectors in terms of complex 
hyperspherical coordinates. Although it is not difficult to find constructive control 
schemes for state transfer based on Lie group decompositions, most schemes do not 
give explicit expressions for the rotation angles ("generalized Euler angles") in the 
factorization, and thus the rotation angles usually have to computed numerically. By 
parametrizing the initial and target states in terms of hyperspherical coordinates, we 
obtain a factorization where all generalized Euler angles are given explicitly in terms 
of the hyperspherical coordinates of the initial and target states, eliminating the need 
for numerical calculation of the generalized Euler angles, aside from computation of 
the hyperspherical coordinates, which is trivial in terms of computational overhead. 

The factorization is applicable given controls capable of implementing phase 
rotations and population rotations (of cither X or Y type) on a collection of two- 
dimensional subspaces, similar to the general requirements for constructive geometric 
control schemes. Compared to control schemes based on the standard factorization, 
this scheme has the additional advantages that all initial and final phase rotations can 
be combined in a single step and executed concurrently, reducing the time required to 
achieve the state transfer. As with all bang-bang control schemes based on Lie group 
decompositions, the factorization only determines the sequence in which the controls 
are applied and the pulse area (rotation angle) of the control pulses, leaving us with 
considerable freedom to choose the pulse shapes and amplitudes, which can be used 
to further optimize a performance index. Here we have considered optimization of the 
pulse amplitudes for piecewise constant controls such as to minimize a time-energy 
performance index that takes into account the competing goals of trying to minimize 
the transfer time and energy cost of the controls. 

The scheme can be generalized to realize unitary operators. By expressing the 
eigenvectors of the target gate U in hyperspherical coordianates we obtain an explicit 
decomposition for arbitrary unitary operators. Aside from giving explicit expressions 
for the Euler angles in terms of hyperspherical coordinates an advantage of the 
decomposition is that it separates the elementary rotations in such a way as to allow 
concurrent implementation of subsets of operations, which can reduce the control time. 
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Specifially, unlike in many standard decomposition schemes, the Z n and Y n rotations 
do not occur in an alternating sequence but are clustered. Since the Z n rotations are 
mutually commuting, this allows concurrent implementation of many operations. 

Acknowledgment s 

This work is support by National Natural Science Foundationof China (Grant No. 
60974037 & Grant No. 11074307). SGS acknowledges funding from EPSRC ARF 
Grant EP/D07192X/1 and Hitachi. 

References 

[1] H. Rabitz, R. de Vivie-Riedle, M. Motzkus, K. Kompa, Whither the Future of Controlling 

Quantum Phenomena? Science 288 (2000) 824-828. 
[2] M. A. Nielsen, I. L. Chuang, Quantum Computation and Quantum Information, Cambridge, 

Cambridge University Press 2000. 
[3] G. M. Huang, T. J. Tarn, J. W. Clark, On the Controllability of Quantum Mechanical Systems, 

J. Math. Phys. 24 (1983) 2608-2618,. 
[4] C. K. Ong, G. Huang, T. J. Tarn, J. W. Clark, Invertibility of Quantum Mechanical Control 

Systems, Math. Sys. Theor. 17 (1984) 335-350. 
[5] J. W. Clark, C. Ong, T. J. Tarn, G. M. Huang, Quantum non-Demolition Filters, Math. Sys. 

Theor. 18 (1985) 33-53. 

[6] A. Blaquiere, S. Diner, G. Lochak, (edit) Information Complexity and Control in Quantum 

Physics, Springer- Vcrlag, New York, 1987. 
[7] A. P. Peirce, M. A. Dahleh, H. Rabitz, Optimal control of quantum-mechanical systems: 

Existence, numerical approximation, and applications, Phys. Rev. A 37 (1988) 4950-4964. 
[8] D. D'Alessandro, Introduction to quantum control and dynamics, CRC Press, 2007. 
[9] H. M. Wiseman, G. J. Milburn, Quantum Measurement and Control, Cambridge University 

Press, Cambridge, 2010. 

[10] J. Janszky, P. Domokos, S. Szabo, P. Adam, Quantum-state engineering via discrete coherent- 
state superpositions, Phys. Rev. A 51 (1995) 4191-4193. 
[11] Y. Makhlin, G. Schon, A. Shnirman, Quantum-state engineering with Joscphson-junction 

devices, Rev. Mod. Phys. 73 (2001) 357-400. 
[12] C. Rangan, A. M. Bloch, C.Monroe, P. H. Bucksbaum, Control of Trapped-Ion Quantum States 

with Optical Pulses, Phys. Rev. Lett. 92 (2004) 113004. 
[13] D. C. Brody and D. W. Hook, On optimum Hamiltonians for state transformations J. Phys. A: 

Math. Gen. 39 167-170 (2006); ibid. 2007, 40 10949. 
[14] S.Y. Lee, H. Nha, Quantum state engineering by a coherent superposition of photon subtraction 

and addition, Phys. Rev. A 82 (2010) 053812. 
[15] A. C. Doherty, K. Jacobs, Feedback control of quantum systems using continuous state 

estimation, Phys. Rev. A 60 (1999) 2700-2711. 
[16] R. van Handel, J. K. Stockton, H. Mabuchi, Feedback control of quantum state reduction, IEEE 

Trans. Automat. Contr. 50 (2005) 768-780. 
[17] N. Yamamoto, K. Tsumura, S. Hara , Feedback control of quantum entanglement in a two-spin 

system, Automatica 43 (2007) 981-992. 
[18] S. G. Schirmer, X. Wang, Stabilizing open quantum system by Markovian reservoir engineering, 

Phys. Rev. A 81 (2010) 062306 
[19] F. Ticozzi, S. G. Schirmer, X. Wang, Stabilizing generic quantum states with Markovian 

Dynamical Semigroups, IEEE Trans. Autom. Control 55(12) (2010) 2901-2905 
[20] J. Zhang, R.-B. Wu, C.-W. Li, T.-J. Tarn, Protecting coherence and entanglement by quantum 

feedback controls, IEEE Trans. Automat. Contr. 53 (2010) 619-633. 
[21] X. Wang, A. Bayat, S. Bosc and S. G. Schirmer, Global Control Methods for GHZ State 

Generation on a ID Ising Chain, Phys. Rev. A 82 (2010) 012330. 
[22] X. Wang, A. Bayat, S. G. Schirmer, S. Bosc, Robust Entanglement in Anti-fcrromagnctic 

Heisenberg Chains by Single-spin Optimal Control, Phys. Rev. A81 (2010) 032312 
[23] M. Reck, A. Zeilinger, H. J. Bernstein, P. Bcrtani, Experimental realization of any discrete 

unitary operator, Phys. Rev. Lett. 73 (1994) 58-61. 



Bang-Bang Control Design for Quantum State Transfer 



17 



[24] V. Ramakrishna, K.L. Flores, H. Rabitz, R. J. Ober, Control of a coupled two-spin system 

without hard pulses, Phys. Rev. A 62 (2000) 053409. 
[25] D. D'Alessandro, The optimal control problem on SO(4) and its applications to quantum control, 

IEEE Trans. Automat. Contr. 47 (2002) 87-92. 
[26] S. G. Schirmcr, A. D. Greentree, V. Ramakrishna, and H. Rabitz, Constructive control of 

quantum systems using factorization of unitary operators, J. Phys. A 35 (2002) 8315-8339. 
[27] R. Cabrera, T. Strohecker, and H. Rabitz, The canonical coset decomposition of unitary matrices 

through Householder transformations, J. Math. Phys. 51 (2010) 082101. 
[28] R. Romano, D. D'Alessandro, Environment-Mediated Control of a Quantum System, Phys. Rev. 

Lett. 97 (2006) 080402. 

[29] M. Zhang, H. Y. Dai, X. C. Zhu, X. W. Li, and D. Hu, Control of the quantum open system by 

quantum generalized measurement Phys. Rev. A 73 (2006) 032101. 
[30] X. Wang, S. G. Schirmer, Arbitrarily-high steady-state entanglement between non-interacting 

atoms via collective decay and adiabatic control (2010), arXiv:1005.2114 
[31] H. M. Wiseman, G. J. Milburn, Squeezing via feedback, Phys. Rev. A 49 (1994) 1350-1366. 
[32] M. Yanagisawa, H. Kimura, Transfer Function Approach to Quantum Control-Part I: Dynamics 

of Quantum Feedback Systems, IEEE Trans. Automat. Contr. 48 (2003) 2107-2120. 
[33] M. Mirrahimi, P. Rouchon, G. Turinici, Lyapunov control of bilinear Schrodinger equations, 

Automatica 41 (2005) 1987-1994. 
[34] M. R. Jovanovic and B. Bamieh, Lyapunov-based distributed control of systems on lattices, 

IEEE Trans. Automat. Contr. 50 (2005) 422-433. 
[35] C. Altafini, Feedback stabilization of isospectral control systems on complex flag manifolds: 

Application to quantum ensembles, IEEE Trans. Automat. Contr. 52 (2007) 2019-2028. 
[36] S. Kuang, S. Cong, Lyapunov control methods of closed quantum systems, Automatica 44 (2008) 

98-108. 

[37] X. Wang, S. G. Schirmer, Entanglement generations between distant atoms by Lyapunov 

Control, Phys. Rev. A 80 (2009) 042305. 
[38] X. Wang, S. G. Schirmer, Analysis of Lyapunov method for control of quantum states, IEEE 

Trans. Atom. Control 55 (2010) 2259-2270; Analysis of Effectiveness of Lyapunov Control for 

Non-generic Quantum States, IEEE Trans. Autom. Control 55 (2010), 1406-1411. 
[39] Shai Machnes, U. Sander, S. J. Glaser, P. de Fouquieres, S. G. Schirmer, T. Schultc-Herbruggen, 

A Unified Programming Platform for Optimising and Benchmarking Quantum Control 

Algorithms (2010) |arXiv:1011.4874| 
[40] D. J. Maas, D. I. Duncan, R. B. Vrijen, W. J. van der Zande, and L. D. Noordam, Vibrational 

ladder climbing in NO by (sub)picosecond frequency-chirped infrared laser pulses, Chem. 

Phys. Lett. 290 (1998) 7580. 
[41] N. Khaneja, R.Brockett, S. J. Glaser, Time optimal control in spin systems, Phys. Rev. A 63 

(2001) 032308. 

[42] L. M. K. Vandersypen and I. L. Chuang, NMR Techniques for quantum control and computation, 

Rev. Mod. Phys. 76 (2004) 1037 
[43] N. Motee, A. Jadbabaie, Optimal Control of Spatially Distributed Systems, IEEE Trans. 

Automat. Contr. 53 (2008) 1616-1629. 
[44] F. Ticozzi, A. Ferrante, M. Pavon, Robust steering of n-level quantum systems, IEEE Trans. 

Automat. Contr. 49 (2004) 1742-1745. 
[45] S. G. Schirmer, Implementation of quantum gates via optimal control, J. Mod. Opt. 56 (2009), 

831. 

[46] R. Nigmatullin, S. G. Schirmer, Implementation of Logic Gates on Encoded Qubits via Optimal 

Control, New J. Phys. 11 (2009) 105032. 
[47] S. G. Schirmer and P. J. Pemberton-Ross, Fast high-fidelity information transfer in spin chain 

quantum wires, Phys. Rev. A 80, 030301 (2009) 
[48] D. D'Alessandro, Optimal evaluation of generalized Euler angles with applications to control, 

Automatica 40 (2004) 1997-2002. 
[49] K. Ch. Chatzisavvas, C. Daskaloyannis, C. P. Panos and S. G. Schirmer, Explicit algorithm for 

generalized Euler angle decomposition of SU(2), Phys. Rev. A 80 (2009), 052329. 
[50] I. Bengtsson, K. Zyczkowski, Geometry of Quantum States: An Introduction to Quantum 

Entanglement, Cambridge, Cambridge University Press, 2006. 
[51] M. A. Nielsen, M. R. Dowling, M. Gu, A. C. Doherty, Quantum Computation as Geometry, 

Science 311 (2006) 1133. 

[52] M. Zhang, H. Y. Dai, H. W. Xie, D. Hu, Controllability of multiple qubit systems, Eur. Phys. 
J. D 45 (2007) 331-334. 



